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Abstract 


This paper presents a comparison between the dynamic behavior of a 250kW stand-alone proton exchange membrane fuel cell power plant 
(PEM FCPP) and a 250kW stand-alone microturbine (MT). Dynamic models for the two are introduced. To control the voltage and the power 
output of the PEM FCPP, voltage and power control loops are added to the model. For the MT, voltage, speed, and power control are used. Dynamic 
models are used to determine the response of the PEM FCPP and MT to a load step change. Simulation results indicate that the response of the MT 
to reach a steady state is about twice as fast as the PEM FCPP. For stand-alone operation of a PEM FCPP, a set of batteries or ultracapacitors is 
needed in order to satisfy the power mismatch during transient periods. Software simulation results are obtained by using MATLAB®, Simulink®, 


and SimPowerSystems®. 
© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


Distributed generation (DG) is one of the most promising 
alternatives for generating and delivering electric power. The 
need for DG is heightened due to the restructuring of the elec- 
tric power industry and the increase in electric power demand. 
In addition, the need becomes more acute because of stringent 
present-day power quality and system reliability requirements. 
Generally, a DG system consists of small-scale power genera- 
tors that are located close to load centers. Primary advantages of 
the DG system are: (a) consumers can generate electric power 
with or without grid backup; (b) the excess generation can be 
sold back to the grid at low load conditions. Fuel cells and 
microturbines as DG sources are among the potential candi- 
dates for supplying electric and thermal energy to residential 
and commercial loads [1]. 

Fuel cells are not only characterized by higher efficiency than 
conventional power plants, but they are also environmentally 
clean, have extremely low emission of oxides of nitrogen and 
sulfur and have very low noise. The main components of the fuel 
cell system include the fuel processing unit or the reformer, the 
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fuel cell stack and the power-conditioning unit [2,3]. Reformer 
produces hydrogen to supply the fuel cell stack by processing 
any hydrocarbon fuel such as propane, methane, or methanol. 
Using the hydrogen input and oxygen from air, fuel cell stack 
produces electricity and water through an electrochemical pro- 
cess. The output from a fuel cell is dc power. To provide power to 
a residential load, or to the electrical grid, a power-conditioning 
unit is needed. Many models have been proposed to simulate the 
fuel cell in the literature [4-8]. Due to the low working temper- 
ature (80-100 °C) and fast start up, proton exchange membrane 
(PEM) fuel cell power plants (FCPPs) are the best candidates 
for residential and commercial applications. Based on the model 
introduced in [4,5], a model for a 250 kW PEM FCPP is devel- 
oped and used to study the dynamic behavior in response to a 
step change. 

Microturbines (MT) are small and simple gas turbines that 
have three main components: compressor, combustor, and tur- 
bine. The high-pressure air from the compressor when mixed 
with the injected fuel forms a combustible mixture. The mixture 
is ignited in the combustor to produce hot gas flow, which is used 
to drive the turbine [9]. Two distinct types of MT are identified in 
the literature, the single- and split-shaft MT [10]. Single-shaft 
MT is high speed in nature, where the electric generator and 
turbine are mounted on the same shaft. Split-shaft MT uses a 
gearbox to connect the electric generator to the turbine shaft. 
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Many models are used in the literature to describe the behavior 
of the MT [10-14]. In this paper, the model used in [11] has 
been modified by adding speed and voltage control loops. 

This paper is organized as follows. Section 2 presents the 
modified 250 kW PEM FCPP model. Section 3 presents the MT 
model. Tests and results are discussed in Section 4. Section 5 is 
the conclusion. 


2. PEM fuel cell model 
2.1. Model description 


In [4,5,15] a model of a PEM FCPP is introduced. In this 
paper the model has been modified to simulate a 250 kW PEM 
fuel cell. 

The model is based on simulating the relationship between 
the output voltage and partial pressures of hydrogen, oxygen, and 
water. A detailed model is shown in Fig. 1, which includes the 
PEM FCPP, reformer, and power conditioning unit models [4]. 
The model is verified based on the experimental results obtained 
using the PEM FCPP, which is currently operated in authors’ 
laboratory. All of the model parameters were tuned to make the 
model produce similar response as the real system. 

The 250kW PEM FCPP model parameters are based on 
440 V dc bus voltage, stack current capacity of 94 A, and cell 
voltage of 0.8 V. Based on the above figures, the PEM FCPP 
consists of six parallel stacks, each stack has 550 cells in series. 
Using the indicated number of cells and stacks the 250 kW PEM 
FCPP model parameters are given in Table 1, where, E is the 
Nernst voltage (V), pp, is the hydrogen partial pressure (atm), 
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Fig. 1. PEM FCPP system block diagram. 


Table 1 

PEM FCPP model parameters 

Parameter Value 

Stack temperature 343K 

Faraday’s constant (F) 96484600 C kmol- ! 
Universal gas constant (R) 8314.47 Jkmol7! K7! 
No load voltage (Eo) 0.8 V 

Number of cells per stack, No 550 

Number of stacks, Netack 6 


K, constant = No/(4F) 
Utilization factor (U) 
Hydrogen valve constant (Ky, ) 
Water valve constant (Ky,0) 
Oxygen valve constant (Ko, ) 
Hydrogen time constant (ty, ) 
Water time constant (tH,0) 
Oxygen time constant (To, ) 
Reformer time constant (T1) 
Reformer time constant (T2) 
Reformer PI gain (C1) 
Conversion factor (CV) 
Activation voltage constant (B) 
Activation voltage constant (C) 
Internal resistance (Ri™) 
External line reactance (X) 

PI gain constants C2, C3 
Voltage reference signal (V;) 


Methane reference signal (Qmethref) 
Hydrogen—oxygen flow ratio (rH_0) 


Current delay time constant (Ty) 


1.4251 x 1076 kmol s~! A7! 
0.8 

4.22 x 1075 kmol s~! atm7! 
7.716 x 10% kmol s~! atm7! 
2.11 x 10% kmol s~! atm7! 

3.37 s 

18.418 s 

6.74s 

2s 

2s 

0.25 

2 

0.04777 A7! 

0.0136 V 

0.2778 Q 

0.05 2 

0.1, 10 

1.0 p.u. 

0.000015 kmol s~! 

1.168 

3s 


pmo İs the water partial pressure (atm), po, is the oxygen partial 
pressure (atm), 10 is the input molar flow of oxygen (kmol s~!), 


and diy is the hydrogen input flow (kmol s~!). 


2.2. FCPP power control 


Power control scheme has been developed in [15] that can be 


summarized as follows: 


Vau V, 
Pao = Zeel "s A S sin(ô) 


(1) 


where Pac is the ac power, m the modulation index, Vee dc 
voltage, V; the load voltage, X the external line reactance, and 6 
is phase angle of the ac voltage mVeel}. 

Assuming a lossless inverter: 


Pac = Pac = Veen! (2) 


_ NstackNol 


2FU H 


qH: 

where Pac is the de power, J the stack current, qu, input molar 

flow of hydrogen, Nstack the number of stacks, No the number of 

cells per stack, F Faraday’s constant, and U is utilization factor. 
From (1), (2), and (3): 


2FUX 


sin(6) = ——___—— qH, 
(8) mVsNoNstack 


(4) 
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Assuming a small phase angle sin(ô) = ô, 


_ 2FUX 


= ———— qu (5) 
mV,NoNstack - . 


Eq. (5) describes the relationship between output voltage phase 
angle ô and hydrogen flow qy,. Eqs. (1) and (5) show that the 
active power as a function of the voltage phase angle 6 can be 
controlled using the amount of hydrogen flow. 

Output voltage can be controlled by the modulation index m. 
The modulation index is controlled using a PI controller. The 
input to the PI controller is the error signal (difference between 
ac terminal voltage Vac and reference voltage V;). 


3. Microturbine model 
3.1. Model description 


As mentioned earlier, microturbines are classified in two cat- 
egories, the single-shaft or high-speed turbine and split-shaft or 
low speed turbine. In the single-shaft configuration, the com- 
pressor, turbine, and electric generator are mounted on the same 
shaft. The turbine speed is in the range of 50,000—120,000 rpm. 
The frequency of the produced voltage will in the range of 
1500-4000 Hz. To reduce the frequency to 60 Hz, a cyclo- 
converter is used. In split-shaft microturbine, the electric 
generator is driven through a gearbox. The gearbox is used to 
reduce the speed to 3600 rpm. Assuming a two-pole synchronous 
generator (SG), the frequency will be 60 Hz. In this case, no 
power electronic devices are needed for frequency conversion. 

In this paper, a split-shaft model is used to determine the 
dynamic behavior of a microturbine. In [11], the authors used 
the GAST model without speed control to simulate the split- 
shaft microturbine. In this paper, due to the use of a SG, a 
speed controller is developed and used with the GAST model. 
Using a two-pole SG with a split-shaft microturbine that runs 
at 3600 rpm, eliminates the need for frequency conversion. The 
model details and parameters are shown in Fig. 2 and Table 2, 
respectively. The focus of this paper is on the electro-mechanical 
behavior of the MT at normal operation condition, where, the 
load can change gradually, or suddenly. Fast dynamics such as 
faults, loss of power, and startup and shutdown transients are not 
considered. 
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Fig. 2. MT system block diagram with speed and power control. 


Table 2 

MT model parameters 

Parameter Value 
Rated power (Prated-MT) 250 kW 
Real power reference (Pref) 1.0 p.u. 


Damping of turbine (Diurbine) 0.03 


Fuel system lag time constant (T1) 10.0s 
Fuel system lag time constant (T2) 0.1s 
Load limit time constant (T3) 3.0s 
Load limit (Lmax) 1.2 
Maximum value position (Vmax) 1.2 
Minimum value position (Vmin) —0.1 
Temperature control loop gain (Kr) 1.0 
Power control proportional gain (Kp) 0.1 
Power control integral gain (Kj) 1.0 
Speed control proportional gain (Ks) 1000 
Speed control integral gain (Ki) 12.5 
Speed reference (Wref) 1.0 p.u. 


K, 
(1+ K, Taos) 


Fig. 3. Synchronous machine block diagram with voltage control. 


3.2. Simplified synchronous generator model 


A simplified model of a synchronous machine is given in 
[16]. A modified version of this model is presented in Fig. 3. 
The equations that are used to drive the synchronous machine 
model are explained in Appendix A. To simulate the syn- 
chronous machine, the authors used a predefined model existing 
in SimPowerSystems® of the MATLAB® software [17]. The 
model parameters are as shown in Table 3, where Ky, Kj are the 
PI gains of the excitation voltage controller (V). 


Table 3 

Synchronous machine model parameters 

Parameter Value 
Rated power (Prated-SG) 250 kW 
Rated line to line voltage (Vratea) 660 V 
Frequency (F) 60 Hz 
Inertia constant (H) 0.822 s 
Damping factor (Kp) 33.63 p.u. 
Number of poles (P) 2 
Internal resistance (R) 0.02 p.u. 
Internal reactance (X) 0.3 p.u. 
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3.3. MT speed control 


As shown in Fig. 2, speed control is achieved by comparing 
the rotor speed œw; with a reference speed wef. A PI controller 
is used to control the error signal (@;e¢ — wr). The PI controller 
output is then connected as a low value gate input. 


3.4. MT real power control 


Referring to Fig. 2, to control the mechanical output power 
from the microturbine, a PI controller is used. The input to the 
PI controller is the error signal (Prep — AP), where AP is the 
difference between the generated and the load power. 


3.5. MT voltage control 


As stated in Table 3, the SG output voltage is 660 V. To 
compare the SG voltage and the PEM FCPP voltage (440 V) 
a transformer is used to step down the SG voltage to 440 V at 
the load terminals. A PI controller is used to control the output 
voltage by controlling the excitation voltage of the synchronous 
generator. The input signal to the PI controller is the differ- 
ence between the output voltage and reference voltage Vier. The 
parameters for the PI controller are Ky =0.005, and K; = 0.1. 


4. Test and results 


To test the fuel cell and the microturbine dynamic models, 
a step change in the load is used as illustrated in Fig. 4. The 
simulation time is 60s, the initial load is 80 kW, and increased 
to 160 kW after 30 s. The effect of stepping the load on the output 
power of the MT and PEM FCPP are shown in Figs. 5 and 6. 

Comparing Figs. 5 and 6, shows that the PEM FCPP took 
about 10s to change the output power to match the load level. 
The delay is mainly caused by the reformer, due to the slow gas 
processing response. On the other hand, the MT has a much faster 
response compared to the PEM FCPP. As shown in Fig. 5, the MT 
response was almost instantaneous, but due to inertia and damp- 
ing factor of MT and SG, the output power oscillates for about 5 s 
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Fig. 4. Step change in load. 
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Fig. 5. MT output power. 
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Fig. 6. PEM FCPP output power. 


before it reaches steady state condition. The above comparison 
showed that the PEM FCPP has poor load following characteris- 
tics. Therefore, for stand-alone operation, the PEM FCPP must 
incorporate a battery or ultracapacitor to satisfy the load require- 
ment during the transient period. Figs. 7 and 8 demonstrate the 
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Fig. 7. Hydrogen flow in PEM FCPP. 
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Fig. 8. Inverter phase angle. 


change in hydrogen flow from the reformer to the stack of the 
PEM FCPP and the phase angle of the inverter. 

Despite the fact that the MT response is much faster than 
that of the PEM FCPP, the mechanical power output from the 
turbine suffers a time delay of about 5s as shown in Fig. 9. 
During this period, the power mismatch is compensated by the 
rotor inertia, which causes a momentarily reduction in the rotor 
speed as illustrated in Fig. 10. During the low rotor speed period, 
the speed controller and the power controller increase the input 
power to the turbine. When the mechanical input power matches 
the electrical power output, the speed controller brings the rotor 
speed back to synchronous speed (Fig. 10). 

The change in the MT and the PEM FCPP voltage at the load 
terminals is demonstrated in Figs. 11 and 12. In Fig. 11, the 
MT voltage is slightly decreased when the load increases due to 
voltage drop in the line and transformer. At t= 30 s, a very short 
period voltage dip of 20 V occurs due to the sudden change of 
the load. As a result of the drop in voltage, the voltage controller 
increases the excitation level by 2.5% to keep the voltage at the 
same level of 440 V (Fig. 13). 
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Fig. 9. Turbine mechanical power. 
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Fig. 12. PEM FCPP voltage at load terminals. 
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Fig. 13. SG excitation voltage. 
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Fig. 14. PEM FCPP stack output voltage. 
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Fig. 15. Inverter modulation index. 


Referring to Fig. 12, the change in the PEM FCPP ac voltage 
is insignificant. This is mainly due to the effectiveness of the 
PEM FCPP voltage controller (Fig. 1). Fig. 14 shows that the 
fuel cell stack output voltage of the PEM FCPP decreases by 
11.82% as the load doubles. The fuel cell stack output voltage 
is fed to a dc/dc boost converter having a constant gain 1.5. The 
change in the modulation index to increase the output voltage 
is inversely proportional to the voltage decrease as illustrated in 
Fig. 15. 


5. Conclusion 


Dynamic behaviors of the PEM FCPP and MT are studied 
and evaluated. Dynamic models of the 250 kW PEM FCPP and 
MT with speed and power control are developed. The step load 
change test results indicate that the response of the MT to reach 
a steady state is about twice as fast as that of the PEM FCPP. 
Due to poor load following characteristics of the PEM FCPP, a 
set of batteries or ultracapacitors are essential in order to satisfy 
the power mismatch during transient period. Although the MT 
response has a better load following capability, the mechanical 
power output experiences a momentary delay. During this delay 
period, rotor speed decreases as part of the rotor kinetic energy 
is used to compensate the power mismatch. By the time the MT 
system reaches a new equilibrium point, the speed control brings 
the rotor speed back to synchronous speed. 

In order to maintain the voltage level of the SG, the voltage 
controller increases the excitation level by 2.5%. Meanwhile, in 
the PEM FCPP, the voltage dip is very insignificant. Although 
its dc output voltage decreases by 11.82%, the change of mod- 
ulation index maintains a constant ac voltage level. 


Appendix A 


Synchronous machine model as explained in [16] is con- 
structed based on synchronous machine swing equation. Swing 
equation for power system dynamics can be expressed in many 
forms. Given % in electrical deg s~!, swing equation is expressed 
as: 


(A.1) 


Linearizing the swing Eq. (A.1) results in an equation as shown 
in (A.2). 
tj Aws = ATm — AT (A.2) 


The basic equations for the simplified linear synchronous gen- 
erator model consist of three equations, i.e., 


K3 K3K4 


AE = A.3 

1 14+ K3Thos 1+ K3Thos ve) 
AT, = K Ad + K? AE, (A.4) 
AV, = K5 AS + Ko AE,. (A.5) 


The constants K1, K2, K3, Ka, Ks, and K6 are dependent upon 
the network parameters, the quiescent operating conditions, and 
the infinite bus voltage. 
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The simplified synchronous generation model as presented in 
Fig. 3 is obtained by combining the effect of the damping torque 
with voltage control, and Eqs. (A.2)-(A.5). The variables are 
defined as follows: 


Tm mechanical torque 

Ta accelerating torque 

Te electromagnetic torque 

Tj time constant (t; = 2H/a,) 

Tao d-axis transient open circuit time constant 

ô rotor angle 

Vi terminal voltage 

Eg rms of the peak stator voltage 

Erp stator EMF 

Or rotor angular speed 

Kı change in the electrical torque for small change in rotor 
angle at constant d-axis flux linkage 

K2 change in the electrical torque for small change in the 
d-axis flux linkage at constant rotor angle 

K3 impedance factor 

K4 constant related to demagnetizing effect 

K5 change in the terminal voltage V; for a small change in 
rotor angle at constant d-axis flux linkage 

K6 change in the terminal voltage V; for a small change in 
the d-axis flux linkage at constant rotor angle 

Kp damping constant 

H inertia constant. 
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